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The isospin character, the collective or single-particle nature, and the sensitivity to the slope of 
the nuclear symmetry energy of the low-energy isovector dipole response (known as pygmy dipole 
resonance) are nowadays under debate. In the present work we study, within the fully self-consistent 
non-relativistic mean field (MF) approach based on Skyrme Hartree-Fock plus Random Phase Ap- 
proximation (RPA), the measured even-even nuclei ''^^Ni, "^Sn and ^°*Pb. To analyze the model 
dependence in the predictions of the pygmy dipole strength, we employ three different Skyrme pa- 
rameter sets. We find that both the isoscalar and the isovector dipole responses of all three nuclei 
show a low-energy peak that increases in magnitude, and is shifted to larger excitation energies, 
with increasing values of the slope of the symmetry energy at saturation. We highlight the fact 
that the collectivity associated with the RPA state(s) contributing to this peak is different in the 
isoscalar and isovector case, or in other words it depends on the external probe. While the response 
of these RPA states to an isovector operator does not show a clear collective nature, the response to 
an isoscalar operator is recognizably collective, for all analyzed nuclei and all studied interactions. 

PACS numbers: 21.60.Jz, 24.30. Gd, 21.10.Pc, 21.60.Ev, 21.10.Rc, 21.65.Ef, 21.10.Gv, 25.20.-x, 25.60.-t 



I. INTRODUCTION 



Collective phenomena in atomic nuclei have been in 
the past [U, 0] and constitute in the present Q one of the 
most active and interesting topics of research in nuclear 
physics. Experimental data on giant resonances have al- 
lowed us to determine fundamental properties associated 
with the nuclear interaction in the nuclear medium, such 
as the nuclear incompressibility, the isoscalar effective 
mass at saturation, and the nuclear symmetry energy 
at some subsaturation density 0-13 • With the advent of 
new experimental facilities employing rare isotope beams 
(RIBs) p^| - [l^ , the possibility of studying exotic modes 
in unstable nuclei is nowadays feasible. The low-energy 
peak present in the isovector dipole response of proton- 
deficient and neutron-rich nuclei, the so called Pygmy 
Dipole Resonance (PDR) [l3|, has been experimentally 
inv estig ated in several cases such as i7-22q 44,48Qg^ 

llllf , 68Ni [13, 116-120,124,130,132gjj PfTjl^nd 208pb 

20| . In addition, recent extensive theoretical calculations 
indicate that such a low-en ergy peak is a common prop- 
erty of neutron- rich nuclei [2l|. However, for this PDR, 
the isoscalar character, the collective or single particle na- 
ture and the sensitivity to the density dependence of the 
symmetry energy are nowadays under debate [Tsl . [2l| - [27j 
Hereinafter, we will refer to the low-energy peak present 
in the isovector dipole response as pygmy dipole strength 
(PDS) instead of the commonly used PDR since the res- 
onant (collective) nature of such a peak has not been 
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confirmed yet. 

Such an observable is not only important for nuclear 
structure applications. It also impacts on the determina- 
tion of reaction rates in the astrophysical r-process 0]. 
While some theoretical investigations consider the PDS 
as a collective phenomena [Tsl . [2T| - [23j , others ended up 
with opposite conclusions |24| . Within the formers, the 
PDS is basically understood as a resonant oscillation of 
the neutron skin against the isospin saturated proton- 
neutron core. In agreement with this picture, in Ref. (22j 
it has been found that a linear relation between the en- 
ergy weighted sum rule (EWSR) exhausted by the PDS 
and the density dependence of the nuclear symmetry en- 
ergy exists. The symmetry energy is a basic property of 
the nuclear equation of state that plays a crucial role in 
a variety of physical systems. From the very big: the 
size of a neutron star or its composition and structure 
[28l [29| : to the very small: the neutron skin thickness 
of a heavy nucleus |30l433 |. Within the latters, opposite 
to such a picture, the PDS has been postulated just as 
a shell effect dependent on the particular nucleus under 
study: in particular, the authors of Ref. 0] argue that 
the strong fragmentation shown by the strength function 
within their theoretical calculations is indicative of such 
a single particle shell effect [3^. As an example of the 
complexity of the problem, in Ref. [2^ the case of ''^Ca 
have been studied in detail. The authors concluded that 
only one of the low-lying excitations within the energy 
range 5-10 MeV can be described as a pygmy resonance. 
Finally, it is also important to mention that within some 
of the works in which the PDS is assumed to be collec- 
tive, it has been stated that such a low-energy peak in 
the dipole response may, indeed, correspond to a toroidal 
mode 0, 113, HI] or that it is weakly dependent on the 
isovector part of the nuclear effective interaction [2^. 

For our study of the PDS on the measured even-even 
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neutron-rich nuclei ^®Ni, -^^^Sn and ^°^Pb that we take 
as representative of different mass regions, we adopt the 
fully self-consistent non-relativistic mean field (MF) ap- 
proach based on Skyrme Hartree-Fock (HF) plus RPA. 
The MF approach provide a unique framework for the 
study of all nuclei along the periodic table except the 
lightest ones. Such models typically display a rather 
small root mean square (rms) deviation on binding en- 
ergies when compared with a large set of experimental 
data 0, [H, [13 and are able, through the RPA approach, 
to predict the main features of Giant Resonances [H^. 
To assess the sensitivity of our analysis on the nuclear 
model, we employthree Skyrnie parameter sets, namely 
SGII [ail, SLy5 (ai and SkI3 H^. Since the low-energy 
isovector dipole response of neutron-rich nuclei may be 
related with the density derivative of the symmetry at 
saturation, the set of chosen models have been selected 
due to the wide range displayed for their predicted val- 
ues of the L parameter. Such a parameter is defined as 
L = 3poo[9csym(/o)/9p]p^ where Csym(p) is the symmetry 
energy, p is the nucleon density and poo is the nuclear 
saturation density. All the studied nuclei are spherical 
and double-magic. This renders our HF calculations rel- 
atively simple and the analysis clearer since neither pair- 
ing nor deformation should be included. 

First and foremost, we are interested in the theoretical 
study of the main features displayed by the low-energy 
RPA state that give rise to the largest contribution to 
the PDS or. hereinafter, RPA-pygmy state. In our work, 
we shall investigate the isoscalar or isovector character 
displayed by the transition densities associated to the 
RPA-pygmy state, and the most relevant particle-hole 
(jjh) excitations contributing to such a state. In partic- 
ular, we will emphasize that different operators will pro- 
duce a different number of coherent contributions from 
ph amplitudes. This means that in the case of different 
experimental probes one will see the same RPA-pygmy 
state with a different associated degree of collectivity. In 
the final stage of the preparation of this manuscript we 
have become aware that a similar analysis has been per- 
formed in Ref. 

A brief summary of the employed formalism is given in 
Section [ll] where some properties of the interactions we 
use arc also detailed. In Section lllll results are presented, 
analyzed and compared with available experimental data. 
Finally, our conclusions arc laid in the last section. 



II. FORMALISM 

In this section we present the general expression of 
the Skyrme interaction as well as some basic properties 
of the parametrizations used in our analysis. A brief 
description of the RPA formalism is also presented. We 
address the reader to Refs. [s^, [4l|,[43 for further details 
on the Skyrme interaction. 



A. Skyrme interaction 

The Skyrme interaction is a zero-range, velocity- 
dependent interaction that describe nucleons with space, 
spin and isospin variables Yi^cTi and r^. It is commonly 
written as in Ref. [IH, 

y(ri,r2) = t^il + XqP„)5{v) 

+ iti(l + xiP,)[P'25(r) + <5(r)p2] 

+ t2{l + X2Pa)'P' ■ 5{V)-P 
+ it3(l+.T3PaK(RM(r) 

+ iWo(cri +0-2) • [P' X ,5(r)P] , 



(1) 



where r = n - ra, R = i(ri + rz), P = ^(Vi - V2), 
P' is the hermitian conjugate of P (acting on the left), 
= 5(1 + f 1 • (T2) is the spin-exchange operator. 
As mentioned already in the Introduction, we employ 
three Skyrme interactions: SGII [ssf . SLy5 [39| and SkI3 
[40j : as many others, they have been accurately cali- 
brated in order to reproduce some bulk properties (the 
binding energies and charge radii) of few selected stable 
nuclei, as well as some empirical nuclear matter proper- 
ties such as the saturation energy and the saturation den- 
sity (and others depending on the specific set). Through- 
out this work, we are mainly interested on the sensitivity 
of the PDS to the density derivative of the symmetry 
energy at saturation [l^]. Since SGII, SLy5 and SkI3 
are characterized by L equal to, respectively, 37.63 MeV, 
48.27 MeV and 100.52 MeV, they span a quite broad 
range (comparable with the one spanned by most of the 
modern and commonly used MF models available in the 
literature [3l|,[33). 



B. Random Phase Approximation 

The discrete RPA method is well-known from text- 
books m, [il] . In our self-consistent approach, we build 

the residual interaction (V^egijuai) proton-proton 
{qq'=pp), neutron- neutron (qq'=nn) and proton-neutron 
{qq'=pn) channels from the Skyrme-HF energy density 

functional, namely V^cgj^uai = ^^^Jup/SpqSpq'. Then we 
solve fully self-consistently the RPA equations by means 
of the matrix formulation like in Refs. [lH, [ig. One 
should note that the continuum is discretized by setting 
the system in a large box. 

For any operator Fjm the (reduced) transition 
strength or probability is given by 



B{EJ,Q^v) = 



\Fjm 



E 

ph 



X 



ph 



ph 



{p\\Fj\\h) 



-(2) 
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where (z^||Fj||C)) is the reduced matrix element of Fjm 
(see, e.g., Ref. The initial state in all studied nu- 

clei, |0), correspond to the RPA ground-state with zero 
total angular momentum and {v) stands for a generic 
RPA excited state. The latter equation is also written 
in an alternative notation that will turn out to be useful 
for our present purposes. That is, each RPA transition 
|0) — >■ excited via Fjm is composed by all consid- 
ered particle-hole (ph) pairs that couple to a total an- 
gular momentum JAl. The relative contribution of each 
ph excitation to the reduced matrix element (z^||i^j||0) is 

accounted by the X^'^ and 11^^^ RPA amplitudes that 
specify a given eigenvector of the RPA secular matrix 
[43l |. For the analysis of the single particle or collective 
character of a given excitation in the response function, 
it is convenient to write the reduced amplitude as follows: 

A,,iEJ,6^ v) = (X^;;) {p\\Fj\\h). (3) 

This is because Eq. ([3]) allows one to determine the co- 
herency (relative sign) and magnitude {\Aph{EJ,Q — > 
1^)1) of all the ph contributions to the reduced transition 
probability. An RPA state is claimed to be a resonant 
excitation if the corresponding reduced amplitude is com- 
posed by several ph excitations similar in magnitude and 
adding coherently. 

The strength function is defined as usual, 

S{E)^Y.\^v\\FjmH{E-E^), (4) 

where Ei, is the eigenenergy associated to the RPA state 
\v). Its moments can be calculated as 

mfc = IdE E'^SiE) =Y,\{y\\F.m\^El. (5) 

Another quantity of interest that characterizes the re- 
lationship of each excited state with the ground state, 
is the transition density. Its integral with a multi- 
pole operator gives the corresponding transition ampli- 
tude of that operator. With the help of the X'^'"^ and 
Y^'^^ amplitudes of a given RPA state \v)^ one can con- 
struct the radial part of its transition density defined by 
bpy[v) = (z^|/5(r)|6) = bp^{r)Y*j^^(r) as follows, 

^kpWjWi^'—^^ ' (6) 

where Ua(r) is the solution to the Skyrme-HF radial 
equations corresponding to the single particle state a. 
Note that the summations in the expression above can 
be done for neutrons or protons separately. This allows 
one to calculate the neutron and proton transition densi- 
ties bpvq (r) (q = n,p) and define accordingly the isoscalar 
(IS) and isovector (IV) transition densities as 

<5p(,is)(r) = dp,,Xr)+Spu,Ar) and (7) 
5p(iv)(r) ^ Sp,,Xr)-Sp.Jr). (8) 



The interest of the transition densities relies on the fact 
that their spatial shape reveal the nature of the excita- 
tions: volume or surface type, isoscalar or isovector, etc. 
Moreover, they can be used as input in calculations of 
inelastic scattering cross sections. More details of our 
implementation of RPA can be found in Ref. [i^ . 

As our theoretical study will be devoted to the low en- 
ergy dipole response in even-even neutron-rich nuclei, we 
define the isoscalar (IS) and isovector (IV) dipole opera- 
tors ( J = 1) used for the different calculations: 

A 

f['^'> =Y.r^Y,M{n), (9) 

4=1 



A 

Fi';:^^ =Y,nYMn)rA^)- (10) 

i=l 

Note that the lowest order term in the IS dipole operator 
coming from the expansion of the Bessel functions (r) 
does not reproduce a physical excitation but a translation 
of the whole system. This is the reason why the IS dipole 
operator in Eq. Q is proportional to the following term 
(r'^) in such an expansion. 

The translational mode should in principle be decou- 
pled from the physical excitations within the RPA. How- 
ever, as in any numerical implementation, the decoupling 
is not perfect. Therefore, part of this state, known as the 
spurious state, overlaps with the physical RPA states. 
There arc different ways to correct this overlap [i^ . 
Our prescription is detailed in Appendix where we 
show how one can subtract the spurious state from the 
neutron and proton transition densities. The reliability 
of our method can be seen in Fig.[l]where we compare the 
strength function — calculated by convoluting the corre- 
sponding reduced transition probability of (Eq. ^ with 
a Lorcnzian of 1 MeV width — for the isoscalar dipole 
response predicted by the SLy5 interaction [s^ for a test 
nucleus (^°*Pb) in three cases: in one case the spurious 
state has not been subtracted (solid curve) , in the second 
case the spurious state has been subtracted by correcting 
the isoscalar dipole operator Eq. (|9]) with the addition of 
a term —rjriYiMiji) where rj ~ 5(r^)/3 like in Ref. (soj 
(dashed line), and finally in the last case the spurious 
state has been subtracted as explained in Appendix |X] 
(dot-dashed line). 

From this figure one clearly sees that the different pre- 
scriptions for correcting the spurious state are completely 
equivalent. The advantage of our method relies on the 
fact that, by construction, we exactly subtract the spuri- 
ous state from the neutron and proton transition densi- 
ties, and these are among the quantities that we discuss 
in detail below. 

In the case of the IV dipole operator, we have to sub- 
tract the dipole motion associated with the displacement 
of the neutron and proton center-of-mass. It can be done 
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FIG. 1: (Color online) Strength function in the case of the 
SLy5 interaction [s^ for a test nucleus (^'^**Pb) as a function 
of the excitation energy for three cases: (i) the spurious state 
has not been subtracted (solid line), (ii) the spurious state has 
been subtracted by correcting the isoscalar dipole operator 
Eq. ^ (dashed line), and (iii) the spurious state has been 
subtracted as explained in the Appendix 1X1 (dot-dashed line). 



by modifying the operator as follows, 



N A 

^5]rpYiM(rp). (11) 



Last but not least, calculations beyond RPA may of 
course to some extent change the quantitative picture. 
It is known that correlations associated with coupling 
with 2 particle- 2 hole (2p-2h) configurations, or particle- 
phonon configurations, tend to shift the strength down- 
ward and increase its fragmentation. However, at least in 
the cases that have been studied in Refs. [5l|, [s^, these 
changes do not destroy the qualitative features associated 
with the PDS, namely the isospin character, the overall 
behavior of the transition densities and the collective or 
single-particle character of the states. We leave aside in 
our study the debated case of Ca isotopes discussed in 
Refs. iilEllsl. 



III. RESULTS 

In this Section we present a detailed study of the dipole 
response of ^^Ni, ^'^^Sn and ^"^Pb as predicted by the 
Skyrme interactions SGll, SLy5 and Skl3 within the for- 
malism described in Section|lIl The results are organized 
as follows. We present for each nucleus the isoscalar and 
isovector dipole strength functions. 

In order to have a simple estimate for the collectiv- 
ity displayed by the different dipole responses, we plot 



also the reduced transition probabilities in single parti- 
cle units (s.p.u., or Weisskopf units H^j). Such unit is 
based on a macroscopic approach. One evaluates the av- 
erage transition rate of a typical excitation in terms of 
the angular momentum carried by the probe, and the 
mass radius of the nucleus under analysis; in this way, 
the result is nucleus- independent. By following Ref. |47| . 
one can calculate the isovector dipole response in Weis- 
skopf units accounting for the center of mass correction 
as. 



^) for protons 



(— ^) for neutrons 



(12) 



where the radius R is taken to be tqA^^^ and rp is the 
radius of the average sphere that one nucleon occupies, 
at the standard saturation density of 0.16 fm^^ (that is. 
To = 1.14 fm). For the isoscalar dipole case, once the 
spurious state has been subtracted as explained in the 
previous section, one finds 



1 

47r V2 



,\r 



iR^ 



(13) 



with 7] = and (r^) ~ f-R^- The sum over all nu- 

cleons in Eqs. (|12p and ([T51) coincides in good approxi- 
mation (around 10%— 20%) with the corresponding total 
RPA strength, although there exist some mass depen- 
dence: heavier nuclei are better reproduced by the esti- 
mate provided by the Weisskopf units. Such a unit al- 
lows us to account qualitatively for the nature of different 
excitations since a given RPA state will contribute with 
several single particle units if it is collective. Moreover, it 
also enables the comparison between the results obtained 
for different nuclei. 

Then, we focus on the low-energy region in order to in- 
vestigate the RPA-pygmy state leading to the PDS. We 
show, first, the neutron and proton transition densities 
associated with such a state for each interaction and nu- 
cleus. The analysis of the transition densities may be 
very illustrative since they allow one to distinguish some 
spatial details related to the dynamics of every excita- 
tion. For example, one could understand if either nucle- 
ons from the surface or from the interior of the nucleus 
are contributing more to the excitation, and this is cru- 
cial to estimate which reaction is more efficient in excit- 
ing this mode. Besides this, the comparison between the 
neutron and proton transition densities informs us about 
the relative motion of neutrons with respect to protons, 
or in other words, on the isoscalar or isovector character 
of each RPA state. 

For this aim, we also use a local criterion to study 
quantitatively the isoscalar and isovector splitting of the 
RPA-pygmy state [s^ based on the following. At each 
radial distance r,; , where i = 1 . . .N at which the neu- 
tron and proton transition densities are calculated, we 
define that a certain RPA state is 70% isoscalar if at 
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least the 70% of the calculated points fulfill the condi- 
tion > |(Sp[/^''(r)|. Moreover, we will exploit 
the possibility of analyzing the isoscalar or isovector na- 
ture of the RPA-pygmy state in different regions of the 
nucleus. Specifically, we impose the above defined crite- 
ria of isoscalarity in two additional regions: one in the 
internal part of the nucleus, i.e. from fm to R/2 and 
the other in external part of the nucleus, namely from 
R/2 to R. 

Finally, we analyze the most relevant particle-hole 
(jjh) excitations contributing to the RPA-pygmy state. 
To this end, we calculate the magnitude and sign (co- 
herency) with which each ph excitation contribute to the 
isovector (IV) and isoscalar (IS) dipole reduced transi- 
tion probability {B{E1;^) where ^ =IV or IS). For that, 
we have used the isoscalar A'^f^{El; ^ = IS) and isovector 
A'if^{El;^ = IV) reduced amplitudes defined in Eq. 



A. Strength functions in ''^^Ni, ^''^Sn and ^°^Ph 

In this subsection we analyze the main features dis- 
played by the strength function Eq. (jU associated to 
the isovector and isoscalar dipole response of ^^Ni, -'^'^^Sn 
and ^°^Pb. It has been calculated by convoluting the 
corresponding reduced transition probability, Eq. 
with a Lorenzian of 1 MeV width. The low-energy 
dipole response of all studied nuclei has been measured 

TABLE I: Excitation energy E and isoscalar (^ — IS) and 
isovector (^ = IV) reduced transition probabilities _B(i51;^) 
corresponding to the RPA-pygmy states of ®*Ni, ^'^^Sn and 
^"^Pb as predicted by SGII, SLy5 and SkI3 interactions. 





force 


E 
[MeV] 


B{E1-IS) B{E1;IV) 
[fm«] [fm^] 


68Ni 


SGII 
SLy5 
Skis 


9.77 
9.30 
10.45 


1.9x10^ 
1.7x10^ 
3.0x10^ 


1.4 
0.8 
3.6 


"^Sn 


SGII 
SLy5 
SkI3 


8.52 
8.64 
9.23 


3.3x10^ 
1.0x10"* 
1.1x10'* 


1.2 
1.6 
7.4 




SGII 
SLy5 
SkI3 


7.61 
7.74 
8.01 


1.7x10" 
2.8x10* 
1.9x10'' 


2.9 
2.8 
6.6 



We start analyzing the results for ^^'^Pb. In Fig. [2^ we 
show the strength function corresponding to the isovec- 
tor dipole response as a function of the excitation energy. 
The inset displays in a larger scale the pygmy region. In 
Fig. ^jp, the same quantities are shown but this time for 
the isoscalar dipole response as a function of the excita- 
tion energy. In both figures, the predictions of the three 




E (MeV) 

FIG. 2: (Color online) Strength function corresponding to the 
isovector (a) and isoscalar (b) dipole response of Pb as a 
function of the excitation energy. The inset in (a) displays in a 
larger scale the pygmy region. In both figures the predictions 
of SGII, SLy5 and SkI3 are depicted. Black arrows indicate 
the experimental centroid energies for the PDS {E = 7.37 
MeV within a window of 6 - 8 MeV) for the ISGDR 

(E = 20.3 ± 2 MeV [Ml) and the energy peak for the IVGDR 
{E = 13.43 MeV and a total width of 2.42 MeV [13]). 



selected interactions are shown. The centroid energies of 
the PDS and the Isoscalar Giant Dipole Resonance (IS- 
GDR) as well as the energy peak in the Isovector Giant 
Dipole Resonance (IVGDR) for ^"'^Pb as predicted by the 
employed interactions {E = 7.6 - 8.0 MeV, £: = 20 - 21 
MeV and = 12 — 13 MeV, respectively) fairly agree 
with the experimental data {E = 7.37 MeV within a 
window of 6 - 8 MeV E = 20.1 - 20.5 MeV [s^ 
and E = 13.43 MeV and a total width of 2.42 MeV 
[131, respectively). Consequently, the RPA predictions 
of SGII, SLy5 and SkI3 may allow us to elucidate the 
microscopical structure and properties of the PDS. In 
Table [Tj the excitation energy and isoscalar and isovector 
reduced transition probabilities of the RPA-pygmy state 
— i.e. the RPA state which give rise to the largest peak 
in the PDS region — are detailed for all the sudied nu- 
clei as predicted by SGII, SLy5 and SkI3. In the case 
of ^°*Pb we find an excitation energy of E = 7.61 MeV 
for SGII, E = 7.74 MeV for SLy5 and E = 8.01 McV 
for SkI3. We qualitatively observe that the low-energy 
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E (MeV) 

FIG. 3: (Color online) Same as Fig. [2] for ^^^Sn. The experi- 
mental value for the peak energy of the PDS {E = 9.8 ± 0.7 
MeV) is indicated by a black arrow [Ts! ]. 



peak found in the IV and IS dipole responses of ^"®Pb 
shows an increasing and outward trend with the excita- 
tion energy as the value of the parameter L increases. 
This behavior is in agreement with Ref. [22j where the 
energy weighted sum rule or mi for the PDS was found 
to be linearly correlated with L in mean-field models. 

In the case of ^^^Sn and ^^Ni, the strength functions 
for the dipole response are depicted in Figs. and 0^ 
(IV) and Figs. ISb and|3}D (IS), respectively. Again, the 
predictions of SGII, SLy5 and SkI3 {E = 8.5 - 9.2 MeV 
for i32Sn and = 9.3 - 10.4 MeV for ^^Ni) are in rather 
good agreement with the measured data {E = 9.1 — 10.5 
MeV for "^Sn ^ E = II MeV and an eiie^gy 
width estimated to be less than 1 MeV for ^^Ni p7|). 
In the case of ^^^Sn, the RPA-pygmy state predicted by 
SGII correspond to a state with an excitation energy of 
E = 8.52 MeV while SLy5 predicts E = 8.64 MeV and 
Skis E = 9.23 MeV. And for ^^Ni, the values predicted 
by SGII, SLy5 and SkI3 for the excitation energy of the 
RPA-pygmy state are respectively: E = 9.77 MeV, E = 
9.30 MeV and E = 10.45 MeV. Qualitatively in both 
nuclei, it seems again that the larger the value of L, the 
higher the values predicted for the excitation energy and 
the larger the different peaks arising in the low-energy 
region (see Figs. |3]and|3]). In addition, we observe for 
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FIG. 4: (Color online) Same as Fig. [2] for "^Ni. The experi- 
mental value for the peak energy of the PDS {E = 11 MeV 
and an energy width estimated to be less than 1 MeV) is 
indicated by a black arrow [ij. 

all nuclei that the PDS is an order of magnitude smaller 
than the IVGDR and that its isoscalar counterpart is 
of the same order of magnitude than the corresponding 
ISGDR. 



B. Reduced transition probability in single particle 
units as an indicator of collectivity 

In Fig. [5l we focus on the relevant region for the study 
of PDS and show the reduced transition probabilities in 
single particle units [see Eqs. ([T^ and and the re- 
lated discussion]. The excitation energies of the RPA- 
pygmy state are also depicted. We display again both the 
isovector (Figs.[S^,[5J: and[Sj3) and isoscalar (Figs. [SJd, [SJi 
and[nf) dipole responses. 

Firstly, we focus on the isovector dipole response of 
208pj-, j-ggg pjg Q^j. calculations predict an RPA- 

pygmy state characterized by « 2-4 single particle units: 
this result does not pin down clearly the nature of the 
state. As a reference, the RPA state leading to the 
largest values of the reduced transition strength in the 
IVGDR contribute with about 30 single particle units if 
the strength is fragmented, and with more than 60 if the 
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for our study of the RPA-pygmy state. 



strength is concentrated in one single peak. This is a 
clear indication of the collective nature of the IVGDR. 
From Fig. [SJd, where the isoscalar or compression dipole 
response of the same nucleus is depicted for the used 
Skyrme interactions, the RPA state leading to the pygmy 
peak is contributing with 15-20 single particle units, very 
similarly in magnitude to those displayed by the largest 
peak in the same isoscalar response at larger excitation 
energies and that can be seen in Fig.l^b. These large val- 
ues indicate the collective character of the RPA-pygmy 
state when it is excited by an isoscalar probe. 

In Figs. Et andEJ; (IV) and Figs. EH andEf (IS), the 
reduced transition probabilities in single particle units 
for the case of ^^^Sn and ^^Ni are depicted, respectively. 
Note that we show only the low energy region. As men- 
tioned, the single particle units normalize the absolute 
value of B{E1;£^) for all nuclei. This statement is qual- 
itatively fullfiled by all the interactions: B{E1, IV) cor- 
responding to the RPA-pygmy state in ^^Ni, ^'^^Sn and 
208pb as calculated with SGII lead to 1.9, 1.1 and 1.9 s.p. 
units, respectively. For the case of SLy5, we find 1.1, 1.5 
and 1.9 s.p. units, respectively. And, finally, for the case 
of SkI3, we find 4.9, 6.4 and 4.2 s.p. units, respectively. 
For the isoscalar response, the B{E1, IS) of the same ex- 
cited states in ''^Ni, "^Sn and ^^^Pb as calculated with 
SGII lead to 12, 9 and 12 s.p. units, respectively. For 
the case of SLy5, we find 11, 18 and 19 s.p. units, re- 
spectively. And, finally, for the case of SkI3, we find 20, 
19 and 14 s.p. units, respectively. 



Despite the fact that the reduced transition probabil- 
ities in s.p. units can be only used as a qualitative es- 
timator of the collectivity displayed by a given RPA ex- 
citation, we can conclude that while the isoscalar dipole 
response of all studied nuclei and within all employed in- 
teractions seems to indicate that the RPA-pygmy state 
develope a certain amount of collectivity, the isovector 
response of the same excited state does not provide a 
clear trend: the collectivity displayed is very weak and 
depends on the used model. 



C. The low-energy RPA states: isoscalar or 
isovector character? 

A collective excited state can be said to be purely 
isovector if the transition densities of protons and neu- 
trons scale as Z and N, respectively, and have opposite 
phase. On the other side, a collective excited state can be 
defined as purely isoscalar if the neutron and proton tran- 
sition densities scale in absolute value in the same way, 
but have the same sign. These two cases are extreme. 
As isospin is not a good quantum number in finite nu- 
clei, and is more and more broken as the neutron excess 
increases, the most common situation corresponds to a 
mixture of a certain degree of isoscalarity and isovectori- 
ality, that can be better seen by looking at the neutron 
and proton transition densities. 

The isoscalar or isovector nature of the low-energy 
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FIG. 6: (Color online) Neutron and proton [(a), (c) and (e)] and isoscalar and isovector [(b), (d) and (f)] transition densities 
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RPA states has been already studied in Ref. [Ill in 
the case of ^'"'Ce. In this work, it was found that the 
low- lying dipole states of ^*°Ce are split into two groups 
depending on their isospin structure. More recently, sim- 
ilar conclusions were found in a study of the pygmy dipole 
strength in ^^^Sn [H^, where it has been stated that the 
theoretical calculations were dominated by a low-lying 
isoscalar component basically due to oscillations of the 
neutron skin thickness of the nucleus under study. It is 
important to note that both investigations were reported 
to be in qualitative agreement with the available experi- 
mental data. 

On the basis of the above mentioned works and the 
definitions given in Eqs. ([5|) and in the text around, we 
present a more systematic study of the isospin struc- 
ture of the low-energy RPA states as predicted by the 
forces SGII, SLy5 and SkI3 for the studied ^SNi, ^^^Sn 
and ^°*Pb nuclei. First of all, we plot the neutron and 
proton, as well as the isoscalar and isovector transition 
densities corresponding to the RPA-pygmy state as pre- 
dicted by each interaction in order to illustrate how these 
low-energy transition densities behave. 

We show the neutron and proton transition densities 
in Figs. ^ and [Bt, and the isoscalar and isovector 
transition densities in Figs. [BJi and |Bf, respectively. 



All of them, correspond to the RPA-pygmy state. The 
position of the proton (rp) and neutron (r„) rms radii 
corresponds to the edges of the grey region that defines 
in this way the neutron skin thickness predicted by each 
interaction. 

For the case of ^°^Pb. it can be seen from Fig. [5^. that 
neutrons and protons oscillate differently depending on 
the interaction but in all cases the surface has a domi- 
nant isoscalar character. On the contrary, the interior or 
bulk region is not dominated by the isoscalar or isovector 
component but it is a mixture of them. The isoscalar or 
isovector dominance is better seen in Fig.|B]3. At the sur- 
face of the nucleus the isovector transition density of the 
RPA-pygmy state is very close to zero, while the isoscalar 
one is not. 

In Figs.|6}; (protons and neutrons) and[6ji (IS and IV), 
we display the transition densities for the case of ^■^^Sn. 
It is interesting to note that the situation is very similar 
to the one found in ^"^^Pb. 

The neutron and proton transition densities corre- 
sponding to the RPA-pygmy state in ^*Ni are depicted 
in Fig. [6^, and the corresponding isoscalar and isovector 
ones are displayed in Fig.[Bf. The behavior of the different 
transition densities is predicted to be very similar within 
the studied models. This did not hold for ^-^^Sn and ^°^Pb 
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where some qualitative differences arose. Therefore, it is 
even more clear in this case that the interior of ^®Ni is 
not dominated by isoscalar or isovector components. At 
the surface of the nucleus, the isoscalar part dominates 
but the isovector part is not very small as it happened 
for i32Sn or ^ospb. 

Then, we apply our criteria for defining a 70% isoscalar 
RPA state [see text after Eq. ^] to all calculated ex- 
cited states and plot their contribution to the isovector 
dipole strength function. We calculate the same quan- 
tity for different regions. First, we apply the criteria to 
those states that are 70% isoscalar in the region between 
and R, where R = roA^^"^ (left panels in Fig. [7]), then 
to those which are 70% isoscalar in the internal part of 
the nucleus, namely between and R/2 (central panels 
in the same figure), and finally to those which are 70% 
isoscalar in the external part of the nucleus between R/2 
and R (right panels of the same figure). Specifically, in 
Fig. [7^, we show the above mentioned calculations for 
^•^^Pb as predicted by SLy5 (dashed line). As a guidance, 
we also show the total isovector dipole stregth function 
(solid line). The results predicted by the other interac- 
tions in the case of ^''^Pb are very similar and we are not 
showing them. From such a figure, it is evident that the 
RPA states which are mostly isoscalar in the whole re- 
gion [0,R] (left panel in Fig. [7^) and in the external part 
of the nucleus [R/2,R] (right panel in Fig. [7^) are essen- 
tially the same ones since both give rise to almost the 
same contributions: most of the PDS and a small contri- 
bution to the rest of the strength function. In the central 
panel, where we represent those RPA states that are 70% 
isoscalar in the internal region of the nucleus [0,R/2], we 
confirm that there is no essential contribution from them 
to the state giving rise to the PDS. 

In Figs.[7}D and[7J;, we show the same quantities but for 
^^^Sn as predicted by SGII and for ^*Ni as predicted by 
SkI3, respectively. Qualitatively, the same results found 
for ^'^®Pb are now found in these figures. Again, we do 
not show the results for the other studied interactions 
since they are very similar and the same conclusions can 
be drawn. 

Our results indicate that one is allowed to qualitatively 
distinguish the PDS from the IVGDR, and state that 
while the latter strength is basically isovector and in- 
volves the motion of mainly internal nucleons, the former 
is more isoscalar than isovector and involves the motion 
of external nucleons, that arc mainly neutrons in a neu- 
tron rich nucleus. 



D. The most relevant particle- hole excitations in 
the low-energy region 
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FIG. 7: (Color online) Strength function corresponding to 
the isovector dipole response of ^"^Pb (a), -^^^Sn (b) and ^*Ni 
(c) as a function of the excitation energy (solid lines), and 
partial contribution due to those states which are at least 
70% isoscalar (dashed line): for each of the panels, the RPA 
states that are at least 70% isoscalar between and R are 
considered in the panel at the left hand side, those which 
have this feature between and R/2 are considered in the 
central panel, and those which are like that between R/2 and 
R are considered in the panel at the right hand side. See text 
for further explanations. 



The RPA states are build as a superposition of particle- 
hole (ph) excitations. A given RPA state shows a collec- 
tive character under the action of an external operator 
if there are many ph excitations providing non negligible 
contributions (each associated with a transition ampli- 



tude and the sum of X!^"^' and ) that add coherently 



in Eq. ([5]). Therefore, within our approach, the collective 
nature of a peak in the response function of a nucleus is 
always hidden in the RPA states. 

For these reasons, our purpose in this subsection is to 
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analyze the contributions of the different ph excitations 
to the RPA-pygmy states depending on the operator used 
to excite the nucleus. In particular, we study the isoscalar 
and isovector dipole responses of the RPA-pygmy states 
since we are most interested in them. Such a study will 
allow us to quantify the collectivity displayed by each 
RPA state depending on the operator, nucleus and inter- 
action used for the theoretical calculations, and detect if 
there are common trends in the microscopic structure. 

First of all, we show in Fig. UK the neutron (black) and 
proton (red) single-particle MF levels close to the Fermi 
energy for ^"^Pb as predicted by SGII (left panel), SLy5 
(right panel) and SkI3 (middle panel). The proton levels 
show a rather similar ordering and spacing for all the 
interactions. On the contrary, the neutron levels differ 
more when different interactions are compared. 

In Figs. IHJd andlHt, also for the case of ^"^Pb, we show 
all the neutron (black) and proton (red) ph contribu- 
tions to the reduced amplitude ylp,j(i?l;^) [see Eq. 
as a function of their excitation energy for the isovec- 
tor and isoscalar dipole responses, respectively. Both re- 
duced amplitudes have been calculated for the case of the 
RPA-pygmy state predicted by the different MF models. 
Notice that not all contributions can be seen from these 
figures since most of them are very small. 

It is evident from Fig. [SJd that the contributions of 
the most relevant ph excitations to the isovector reduced 
amplitude are only a few in number and there is some 
amount of destructive interference. Accordingly, we have 
seen that their total contribution to the isovector reduced 
transition strength in s.p. units do not clearly exceed 
one. Opposite to that, it is also evident from Fig.[Bl: that 
the contributions of the most relevant ph excitations to 
the isoscalar reduced transition amplitude are basically 
coming from neutron transitions, and that most of them 
add coherently. This is, actually, one of the basic features 
to assert that a given RPA state is collective. 

In Table [TI] in Appendix [B] we show the numerical de- 
tails of the ten neutron and proton ph excitations provid- 
ing the largest contributions to the isovector and isoscalar 
74p^(i?l;^). We also indicate in the figures and table 
the involved single particle orbitals corresponding to the 
largest contributions to the reduced amplitudes. In par- 
ticular, we generally find within all the employed models 
that the dynamics of the low-energy isoscalar dipole re- 
sponse of ^°®Pb seems to be governed by the excitations 
of the outermost neutrons, namely those that form the 
neutron skin thickness of this nucleus. From the analysis 
of the ph contributions, we conclude that while the low- 
energy isoscalar dipole response of ^°^Pb arising form the 
RPA-pygmy state can be considered as a collective mode 
in all studied models, the PDS can not. 

Given the relevance of the analysis of the different re- 
duced amplitudes ^^^(i?!;^), we show the same figures 

also for ^^^Sn and ^*Ni where the same features found 
in 208p|^ pg^jj^ |-,g observed. The neutron and proton sin- 
gle particle levels of these nuclei are displayed in Figs. 
[HK and [TUb. respectively. For the case of ^^^Sn, the pro- 
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Fermi level is indicated by a dashed black line. All ph contri- 
butions to the isovector reduced amplitude corresponding to 
the ^"*Pb RPA-pygmy state as a function of the ph excitation 
energy (c). All studied models are shown. Largest neutron p/i 
contributions are also listed in decreasing order from top to 
bottom. Same as (b) but for the isoscalar reduced amplitude 
(c). 



ton single particle levels display a similar spacing and 
ordering for the different models while the neutron ones 
do not. On the contrary, the ordering and spacing of 
the neutron and proton single particle levels in ^*Ni are 
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FIG. 10: (Color online) Same in Fig. [8] for the case of ''^Ni. 



very similar within all studied forces. This difference do 
not affect the qualitative structure of the isovector and 
isoscalar reduced amplitudes for ^-^^Sn (Figs. [Hb andj^t, 
respectively) and ^^Ni (Figs. [TOb and [TOh. respectively). 
As in the case of ^°^Pb. while the isovector dipole re- 
sponses of these nuclei do not show in the corresponding 
reduced amplitudes more than a few relevant neutron and 
proton ph excitations (with some amount of destructive 
interference) , the isoscalar reduced transition amplitudes 
display a more clear collective behavior of the neutron ph 
excitations coming from the transitions of the outermost 



levels (see Tables Hill and ITVl of the Appendix IB] for some 
numerical details). 



IV. CONCLUSIONS 

The collective or single-particle character of the low- 
energy dipole response has been carefully studied in three 
nuclei, representative of different mass regions and neu- 
tron excess. To that end, we have analyzed within the 
fully self-consistent non-relativistic MF Skyrmc Hartree- 
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Fock plus RPA approach, the measured even-even nuclei 
^*Ni, -'^'^■^Sn and ^°®Pb. In order to investigate the sensi- 
tivity of the low-energy isovector and isoscalar dipole re- 
sponses of the studied nuclei on the slope of the symmetry 
energy, we employed three Skyrme interactions: SGII, 
SLy5 and SkI3. They have been selected because they 
cover a wide range for the predicted values of the L pa- 
rameter. We have qualitatively found that the isoscalar 
as well as the isovector dipole responses for all studied 
nuclei show a low-energy peak in the strength function 
that increases in magnitude and is shifted to larger en- 
ergies with increasing values of L. The behavior of the 
isovector dipole response of these nuclei is in agreement 
with recent findings for which a much larger set of rel- 
ativistic and non-relativistic MF interactions were used 

From our analysis of ^^Ni, ^'^^Sn and ^"®Pb, wc have 
also seen that the collectivity associated with the RPA 
states giving rise to the PDS show up depending on the 
nature of the probe used for exciting the nucleus: in par- 
ticular, there is systematically more collectivity in the 
isoscalar than in the isovector transitions. 

To detail this conclusion more, the studied RPA-pygmy 
states consistently display a strong isoscalar character, 
although a non negligible isovector component is always 
observed. Our results do not support a clear collective 
nature of the isovector response. This is opposite to what 
happens for all studied interactions when the same nuclei 
are excited by an isoscalar dipole operator. In this, the 
low-energy peak appearing in the strength function is 
recognizably collective and basically due to the outermost 
neutrons — in other words, to the neutrons forr ning the 
neutron skin known to be present in ^"^Pb [s^ . l59j and 
likely to developed also in ^^Ni and "'^'^^Sn. 

Therefore, we can conclude that the isoscalar dipole 
oscillations of neutron rich nuclei arising from the RPA- 
pygmy states can be understood within the Skyrme-HF 
plus RPA approach as a collective motion of the outer- 
most neutrons in neutron rich nuclei. Such a statement 
should be confirmed by further experimental investiga- 
tions. 



Appendix A: Spurious state 

We subtract the spurious state from the proton (g = p) 
and neutron [q = n) transition densities dpl{r) [where 



ly characterizes the RPA state, cf. Eq. ([6])] by impos- 
ing, first, that the translational operator which is pro- 
portional to the radial coordinate r does not give any 
finite transition amplitude. This means 



(Al) 



where D denotes the corrected RPA state without spuri- 
ous state contamination. As a second condition, we im- 
pose on these new transition densities that the strength 
of the IVGDR is not modified. That is, we write 



drr^r 



By writing 



5pl = 5pl - a' 



dr ' 



(A2) 



(A3) 



where p^^pir) is the proton (q = p) or neutron (q = n) 
density obtained from the self-consistent HF calculations, 
we find the following solution. 



a 



N J drr^ rSp^, 
Z J drr^ rd 



dr 

-2 rSp,, 



A 



J drr^ 



dr 



(A4) 
(A5) 



where Sp„ = 5p2 + Sp^. 



Appendix B: Most relevant ph contributions 

In this Appendix we give the numerical details of 
the most relevant ph contributions to the isoscalar and 
isovector dipole response of ^^Ni, ^'^^Sn and ^"^Pb as pre- 
dicted by SGII, SLy5 and SkI3 studied in Sec.|llll Specif- 
ically, we show for all nuclei and interactions the ten neu- 
tron and proton ph excitations with larger contributions 
to the reduced amplitude A{E1]^) [see Eq. ([3])]. 
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nances: Fundamental High-Frequency Modes of Nuclear 
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Prog. Phys. 70 691 (2007). 
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Phys. A 591, 435 (1995). 
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Bonche Phys. Rev. C 70, 024307 (2004). 
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TABLE II: The ten largest contributions of tlie proton (q — p) and neutron (g = n) ph excitations to the isovector and isoscalar 
reduced amplitudes, Apf^{El; where ^ = IV and IS, respectively, as predicted by SGII, SLy5 and SkI3 models for ^"^Pb. The 
single particle levels involved in the transitions corresponding to each ph excitation are also indicated. 

SGll SkI3 SL^^5 

Isovector 





E = 7.61MeV 




E = 8.01MeV 




E = 7.74MeV 


Aph 








[fm] 






[fm] 






[fm] 


n 


1*13/2 - 


lil5/2 


-2.87 


1*13/2 


lil5/2 


-1.18 


1*13/2 - 


lil5/2 


-2.02 




3p3/2 - 


^ 3^5/2 


1.31 


3pi/2 - 


4Si/2 


1.02 


3pi/2 - 


> 4Si/2 


0.75 




3pi/2 - 


> 3d3/2 


0.81 


l/*9/2 - 


1*11/2 


-0.54 


3pi/2 - 


> 3^3/2 


0.63 




1/19/2 - 


> 1*11/2 


-0.76 


3pi/2 - 


3d3/2 


0.48 


3P3/2 - 


> 3d5/2 


0.54 




2/7/2 — 


> 2^9/2 


0.70 


3P3/2 - 


-S- 3^5/2 


0.27 


1/19/2 - 


> 1*11/2 


-0.53 




2/5/2 - 


> 2g7/2 


0.53 


3P3/2 - 


4Si/2 


0.26 


3P3/2 - 


> 4Si/2 


0.33 




2/5/2 - 


> 3d3/2 


0.28 


2/5/2 - 


3d3/2 


0.12 


2/5/2 - 


> 3d3/2 


0.26 




3P3/2 ~ 


^ 4s, /2 


0.28 


2/5/2 - 


-> 2(77/2 


0.12 


1/19/2 — 


> 2(^7/2 


0.15 




I/I9/2 - 


> 2p7/2 


0.20 


2/7/2 - 


3^5/2 


0.11 


2/7/2 — 


> 3d5/2 


0.12 




3Pl/2 - 


> 4Si/2 


-0.15 


3P3/2 - 


-5> 3^3/2 


0.05 


2/5/2 — 


> 2c?7/2 


0.11 


P 


I/I1I/2 - 


->■ 1*13/2 


3.37 


l/*ll/2 


->■ 1*13/2 


1.52 


l/*ll/2 - 


1*13/2 


2.94 




2(^5/2 — 


> 2/7/2 


-1.30 


lfl7/2 - 


-7 I/19/2 


0.73 


2d5/2 - 


> 2/7/2 


-0.93 




lff7/2 - 


> l/l9/2 


0.85 


3Si/2 - 


-> 3pi/2 


-0.14 


lfl7/2 - 


> I/I9/2 


0.77 




2^3/2 — 


> 2/5/2 


-0.69 


2d5/2 - 


3p3/2 


-0.10 


3Sl/2 - 


>• 3p3/2 


-0.50 




3Sl/2 - 


> 3p3/2 


-0.38 


2rf5/2 - 


2/7/2 


-0.10 


2(^3/2 - 


> 2/5/2 


-0.34 




3Sl/2 - 


> 3pi/2 


-0.19 


2(^3/2 - 


-> 2/5/2 


-0.09 


2d5/2 - 


> 3p3/2 


-0.21 




lff7/2 - 


> 2/5/2 


-0.19 


2d3/2 - 


-7 3pi/2 


-0.07 


2^3/2 - 


> 3pi/2 


-0.21 




2d5/2 - 


> 3p3/2 


-0.15 


3Si/2 - 


-> 3p3/2 


-0.06 


I57/2 - 


> 2/5/2 


-0.18 




2^3/2 — 


> 3pi/2 


-0.13 


157/2 - 


-> 2/5/2 


-0.04 


3Si/2 — 


^ 3pi/2 


-0.18 




2d5/2 - 


> 2/5/2 


-0.06 


l/ill/2 


^ 2.99/2 


-0.02 


20(3/2 - 


> 3p3/2 


-0.06 












Isoscalar 
















[fm^] 






[fm^] 






[fm'l 


n 


l«13/2 - 


ljl5/2 


28.53 


3pi/2 ' 


4Si/2 


40.13 


3pi/2 - 


> 4Si/2 


27.42 




l/l9/2 - 


> 1*11/2 


12.54 


3pi/2 - 


-> 3^3/2 


11.37 


1*13/2 - 


lil5/2 


17.50 




3pi/2 - 


> 3^3/2 


8.86 


1*13/2 ■ 


-> ljl5/2 


9.80 


3pi/2 - 


> 3d3/2 


11.62 




2/5/2 - 


>• 3^3/2 


8.57 


3p3/2 - 


4Si/2 


9.64 


3P3/2 - 


> 4Si/2 


11.16 




3p3/2 - 


> 3^5/2 


8.36 


1/19/2 - 


-S- 1*11/2 


8.55 


2/5/2 - 


!> 3d3/2 


9.10 




3p3/2 - 


^ 4Si/2 


6.79 


2/5/2 - 


-5- 3^3/2 


4.64 


lhg/2 - 


> 1*11/2 


8.46 




1/19/2 - 


> 2^7/2 


6.03 


3P3/2 - 


-i> 3^5/2 


4.44 


3P3/2 - 


> 3d5/2 


6.14 




2/7/2 — 


> 2gg/2 


-5.69 


2/7/2 - 


-7 3d5/2 


3.54 


l/*9/2 - 


> 2^7/2 


4.51 




3pi/2 - 


i> 4Si/2 


-4.18 


l/*9/2 - 


2^7/2 


1.49 


2/7/2 - 


)■ 3d5/2 


3.49 




2/5/2 - 


> 2^7/2 


-2.33 


3P3/2 - 


-i> 3d3/2 


1.22 


2/5/2 — >■ 3d5/2 


1.87 


P 


lftll/2 - 


^> 1*13/2 


22.10 


l/*ll/2 


-S> 1*13/2 


8.79 


l/lll/2 - 


1*13/2 


17.69 




2^5/2 — 


> 2/7/2 


-9.71 


lff7/2 - 


-7 I/I9/2 


7.69 


lp7/2 - 


)■ I/I9/2 


8.22 




lff7/2 l/ig/2 


9.20 


2rf5/2 - 


3p3/2 


0.82 


2d5/2 - 


> 2/7/2 


-6.31 




2^3/2 — 


> 2/5/2 


-4.69 


2d5/2 - 


2/7/2 


-0.65 


lP7/2 - 


> 2/5/2 


2.81 




I57/2 - 


^ 2/5/2 


2.93 


2rf3/2 - 


-7 3pi/2 


0.62 


2(^3/2 - 


> 2/5/2 


-2.07 




2^3/2 — 


> 3pi/2 


1.11 


l/lll/2 


4ii3/2 


0.60 


2(^3/2 - 


> 3pi/2 


1.86 




2d5/2 — 


> 3p3/2 


1.10 


I57/2 - 


-> 2/5/2 


0.59 


2d5/2 - 


> 3p3/2 


1.70 




lftll/2 ^ 


-> 4ii3/2 


0.96 


l/*ll/2 


2^9/2 


0.56 


l/*ll/2 - 


2C/9/2 


1.55 




lftll/2 - 


^> 2^9/2 


0.95 


2rf3/2 - 


2/5/2 


-0.54 


lff7/2 - 


)> 3/19/2 


1.19 




1/5/2 - 


> 2g7/2 


0.92 


lp7/2 - 


-7 4/l9/2 


0.46 


l/*ll/2 - 


4ii3/2 


1.07 



[8] P. G. Reinhard, Nucl. Phys. A 649, 305c (1999). 
[9] L. Trippa, G. Colo, and E. Vigezzi Phys. Rev. C 77, 
061304(R) (2008). 
[10] I. Tanhihata, Prog. Part. Nucl. Phys. 35, 505 (1995). 
[11] H. Geissel, G. Miizenberg, and R. Riisager, Annu. Rev. 

Nucl. Part. Sci. 45, 163 (1995). 
[12] A. Mueller, Prog. Part. Nucl. Phys. 46, 359 (2001). 



[13] A. Klimkiewicz et al, Phys. Rev. C 76, 051603(R) 
(2007). 

[14] A. Leistenschneider et al. Phys. Rev. Lett. 86 5442 
(2001). 

[15] T. Hartmann et al. Phys. Rev. Lett. 85 274 (2000). 
[16] T. Hartmann et al. Phys. Rev. Lett. 93 192501 (2004). 
[17] O. Wieland et. al., Phys. Rev. Lett. 102 092502 (2009). 
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TABLE III: Same as in Table [III for "^Sn. 





SGII 






SkI3 






SLy5 












Isovector 










E = 8.52MeV 




E = 9.23MeV 




E = 8.64MeV 








[fm] 






[fm] 






[fm] 


n I/111/2 - 


lil3/2 


-2.02 


1^*11/2 


1*13/2 


-1.40 


1^*11/2 - 


1*13/2 


-1.59 


2^5/2 - 


> 2/7/2 


1.08 


3Sl/2 - 


3p3/2 


1.12 


3Sl/2 - 


S> 3p3/2 


0.77 


3Sl/2 - 


> 3p3/2 


0.59 


lff7/2 - 


-i> lhg/2 


-0.50 


2(^3/2 - 


> 3pi/2 


0.41 


2^3/2 — 


> 2/5/2 


0.38 


3Sl/2 - 


-i> 3pi/2 


0.41 


2d3/2 - 


> 3p3/2 


0.27 


lS7/2 1^9/2 


-0.33 


2rf3/2 - 


"> 3pi/2 


0.34 


3Sl/2 - 


> 3pi/2 


0.26 


3Sl/2 - 


> 3pi/2 


0.29 


2(^3/2 - 


-> 2/5/2 


0.25 


157/2 - 


> 1/19/2 


-0.23 


2(^3/2 — 


> 3p3/2 


-0.25 


2d5/2 - 


3p3/2 


0.16 


2d5/2 - 


> 2/7/2 


-0.18 


I57/2 - 


^ 2/5/2 


0.23 


2(^3/2 - 


3p3/2 


0.13 


2(^5/2 - 


> 3p3/2 


0.18 


2^3/2 - 


> 3pi/2 


0.22 


1 n , 

^97/2 - 


^75/2 


0.09 


^97/2 — 


^ ^./5/2 


0.15 


lff7/2 - 


^ 2/7/2 


0.07 


2^5/2 - 


2/7/2 


0.07 


2rf3/2 - 


> 2/5/2 


0.15 


V I59/2 - 


> l/lll/2 


1.84 


I59/2 - 


->■ 1^111/2 


1.63 


lff9/2 ~ 


> l/*ll/2 


1.92 


2P3/2 - 


> 2^5/2 


-0.80 


1/5/2 - 


lff7/2 


0.66 


2p3/2 - 


> 2(^5/2 


-0.59 


1/5/2 - 


> lff7/2 


0.27 


2p3/2 - 


2^5/2 


-0.12 


1/5/2 - 


5- lP7/2 


0.37 


2pi/2 - 


> 2d3/2 


—0.23 


2pi/2 - 


2^3/2 


—0.12 


2pi/2 - 


> 3Si/2 


—0.24 


1/5/2 — 


> 2d3/2 


-0.15 


2p3/2 - 


-i> 3Si/2 


-0.10 


2pi/2 - 


> 2d3/2 


-0.20 


2P3/2 - 


> 3Si/2 


-0.08 


2pi/2 - 


3Si/2 


-0.08 


1/5/2 - 


> 2d3/2 


-0.18 


2P3/2 - 


> 2d3/2 


-0.05 


1/5/2 - 


2d3/2 


-0.05 


2p3/2 - 


> 3Si/2 


-0.17 


lff9/2 - 


^ 2/7/2 


-0.04 


1^9/ 2 - 


2/7/2 


-0.04 


1^9/2 - 


^ 2/7/2 


-0.09 


1/7/2 - 


> 2^5/2 


-0.03 


2p3/2 - 


2d3/2 


-0.03 


1/7/2 - 


^ 2rf5/2 


-0.08 


1/5/2 - 


> 2^5/2 


-0.02 


1/7/2 - 


lff7/2 


0.02 


1/5/2 - 


> 2d5/2 


-0.03 










Isoscalar 














[fm^] 






[fm^] 






[fm^] 


n I/111/2 


1*13/2 


12.42 


3Si/2 - 


3p3/2 


24.45 


3Si/2 - 


> 3p3/2 


14.49 


3Si/2 - 


> 3p3/2 


7.69 


2d3/2 - 


3pi/2 


11.25 


2d3/2 - 


> 3pi/2 


13.55 


2(^3/2 - 


> 3p3/2 


-6.21 


3Sl/2 - 


3pi/2 


10.07 


2rf3/2 - 


> 3p3/2 


8.13 


2^3/2 - 


> 3pi/2 


5.87 


lff7/2 - 




5.35 


1^*11/2 ■ 


1*13/2 


7.18 


I57/2 ~ 


> 2/5/2 


5.37 


1^111/2 


1*13/2 


5.29 


3Si/2 - 


> 3pi/2 


5.95 


3Sl/2 - 


>■ 3pi/2 


4.40 


2d5/2 - 


-i> 3p3/2 


4.49 


2(^5/2 - 


> 3p3/2 


4.61 


157/2 - 


> lhg/2 


3.84 


2^3/2 - 


3p3/2 


4.14 


I57/2 - 


> 2/5/2 


3.76 


2^5/2 - 


> 2/7/2 


-2.76 


lff7/2 - 


2/5/2 


2.19 


1^11/2 - 


3gg/2 


2.71 


lff7/2 - 


> 2/7/2 


2.24 


1^11/2 


3^9/2 


1.84 


lff7/2 - 


> 1/1-9/2 


2.51 


2d5/2 - 


> 3p3/2 


1.15 


2(^3/2 - 


2/5/2 


1.81 


I59/2 - 


> 2/7/2 


1.35 


V I59/2 - 


> 1^11/2 


7.59 


lff9/2 - 


-)■ 1^11/2 


5.75 


I59/2 - 


> l/*ll/2 


7.10 


2P3/2 - 


> 2^5/2 


-3.55 


1/5/2 - 


lP7/2 


4.96 


1/5/2 - 


5- lP7/2 


2.81 


1/5/2 - 


> lS7/2 


2.12 


I59/2 - 


2/7/2 


0.64 


2p3/2 - 


> 2d5/2 


-2.26 


1/5/2 - 


> 2d3/2 


1 AO 


1/5/2 - 


2d3/2 


n /I c 
U.4o 


1/5/2 - 


> 2^3/2 


1 

i.DO 


2pi/2 - 


> 2d3/2 


-0.92 


2p3/2 - 


2d5/2 


-0.47 


lff9/2 — 


> 2/7/2 


1.56 


159/2 - 


^ 2/7/2 


0.64 


2pi/2 - 


2^3/2 


-0.41 


2Pl/2 - 


> 3Si/2 


1.02 


159/2 - 


^ 4/lii/2 


0.41 


I59/2 - 


4/lii/2 


0.40 


1/7/2 - 


> 3^9/2 


0.80 


1/5/2 - 


> 4C/7/2 


0.34 


1/5/2 - 


437/2 


0.38 


I59/2 - 


> 4/111/2 


0.78 


1^3/2 - 


> 2/5/2 


0.33 


I59/2 - 


5/lii/2 


0.35 


1/7/2 - 


> 2d5/2 


0.76 


1/7/2 - 


> 3gg/2 


0.28 


2pi/2 - 


~> 3Si/2 


0.33 


14/2 - 


> 2/7/2 


0.69 



[18] P. Adrich et. ai, Phys. Rev. Lett. 95 132501 (2005). 
[19] S. C. Fultz et ai, Phys. Rev. 186, 1255 (1969). 
[20] N. Ryezayeva et ai, Phys. Rev. Lett. 89, 272502 (2002). 
[21] Tsunenori Inakura, Takashi Nakatsukasa, and Kazuhiro 

Yabana, Phys. Rev. C 84, 021302(R) (2011). 
[22] A.Carbone et ai, Phys. Rev. C 81, 041301 (2010). 
[23] J. Piekarewicz, Phys. Rev. C 83 034319 (2011). 
[24] P.-G. Reinhard and W. Nazarewicz, Phys. Rev. C 81, 

051303 (2010). 



[25] D. Gambacurta, M. Grasso, and F. Catara, Phys. Rev. 

C 84 034301 (2011). 
[26] V.Baran, B. Frecus, M. Colonna, and M. Di Toro, 

arXiv:11 11.6504 Vl. 
[27] E. Yiiksel, E. Khan, and K. Bozkurt, Submited to Nucl. 

Phys. A 

[28] J. M. Lattimer and M. Prakash, Astrophys. J. 550, 426 

(2001); Science 304, 536 (2004). 
[29] X. Roca-Maza and J. Piekarewicz, Phys. Rev. C 78, 
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TABLE IV: Same as in Table [III for ^^Ni. 

SGll SkI3 SL^ 

Isovector 





E = 9.77MeV 


Aph 


E = 10.45MeV 


Aph 


E = 9.30MeV 










[fm] 






[fm] 






[fm] 


n 


1/7/2 - 


I59/2 


-1.34 


2pi/2 


24/2 


1.00 


l./'7/2 - 


1^9/2 


-1.01 




2P3/2 - 


24/2 


0.85 


1/7/2 


159/2 


-0.92 


1/5/2 - 


24/2 


0.93 




1/5/2 - 


^ 2d3/2 


0.71 


2p3/2 


^ 24/2 


0.35 


2p3/2 - 


24/2 


0.45 




2pi/2 - 


2d3/2 


0.46 


1/5/2 


2p7/2 


-0.26 


2pi/2 - 


3Si/2 


-0.12 




1/5/2 - 


-> lff7/2 


-0.15 


2p3/2 


-5> 3Si/2 


0.22 


2pi/2 - 


24/2 


0.10 




2p3/2 - 


3Si/2 


0.13 


1/5/2 


24/2 


0.18 


1/5/2 - 


■> 34/2 


0.09 




1/5/2 - 


■> 2^5/2 


-0.12 


1/5/2 


lff7/2 


-0.13 


1/5/2 - 


2p7/2 


-0.08 




2pi/2 - 


-> 3si/2 


-0.04 


1/5/2 


^ 24/2 


0.11 


1/5/2 - 


■> 24/2 


-0.05 




1 f ,„ 

-1-/5/2 


^ •J"3/2 


0.03 


^Pl/2 


^ ■3*1/2 


-0.09 


^P3/2 


^ ■3*1/2 


0.04 




1/7/2 - 


-> 24/2 


— U.Uo 


1/5/2 


34/2 


n no 


1/5/2 - 


lff7/2 


n no 


P 


1/7/2 - 


I59/2 


1.18 


1/7/2 


lff9/2 


1.04 


1/7/2 - 


-5- lP9/2 


0.64 




2Sl/2 - 


2p3/2 


-0.23 


14/2 


^ 1/5/2 


0.31 


14/2 - 


-> 2pi/2 


-0.10 




1^3/2 - 


-> 2pi/2 


-0.14 


2Sl/2 


2p3/2 


-0.12 


2Si/2 - 


2p3/2 


-0.10 




1^3/2 - 


1/5/2 


0.10 


14/2 


2pi/2 


-0.07 


14/2 - 


-> 1/5/2 


0.09 




1/7/2 - 


-> 24/2 


-0.07 


1/7/2 


24/2 


-0.03 


1/7/2 - 


■> 24/2 


-0.06 




1^3/2 - 


2p3/2 


-0.06 


14/2 


2p3/2 


-0.02 


14/2 - 


-> 2p3/2 


-0.03 




2Sl/2 - 


2pi/2 


-0.05 


1^5/2 


^ 1/5/2 


0.01 


2Si/2 - 


2pi/2 


0.03 




1^5/2 - 


2p3/2 


-0.03 


1/7/2 
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-5- 4.99/2 
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